\(\int \frac {(A+B \cos (c+d x)+C \cos ^2(c+d x)) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx\) [1338]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (warning: unable to verify)
   Maple [B] (verified)
   Fricas [A] (verification not implemented)
   Sympy [F(-1)]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 45, antiderivative size = 257 \[ \int \frac {\left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\frac {\sqrt {2} (A-B+C) \arctan \left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}}{\sqrt {a} d}-\frac {2 (43 A-91 B+35 C) \sqrt {\sec (c+d x)} \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d \sqrt {a+a \cos (c+d x)}}+\frac {2 A \sec ^{\frac {7}{2}}(c+d x) \sin (c+d x)}{7 d \sqrt {a+a \cos (c+d x)}} \]

[Out]

2/105*(31*A-7*B+35*C)*sec(d*x+c)^(3/2)*sin(d*x+c)/d/(a+a*cos(d*x+c))^(1/2)-2/35*(A-7*B)*sec(d*x+c)^(5/2)*sin(d
*x+c)/d/(a+a*cos(d*x+c))^(1/2)+2/7*A*sec(d*x+c)^(7/2)*sin(d*x+c)/d/(a+a*cos(d*x+c))^(1/2)+(A-B+C)*arctan(1/2*s
in(d*x+c)*a^(1/2)*2^(1/2)/cos(d*x+c)^(1/2)/(a+a*cos(d*x+c))^(1/2))*2^(1/2)*cos(d*x+c)^(1/2)*sec(d*x+c)^(1/2)/d
/a^(1/2)-2/105*(43*A-91*B+35*C)*sin(d*x+c)*sec(d*x+c)^(1/2)/d/(a+a*cos(d*x+c))^(1/2)

Rubi [A] (verified)

Time = 1.05 (sec) , antiderivative size = 257, normalized size of antiderivative = 1.00, number of steps used = 8, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.133, Rules used = {4306, 3122, 3063, 12, 2861, 211} \[ \int \frac {\left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\frac {\sqrt {2} (A-B+C) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \arctan \left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {\cos (c+d x)} \sqrt {a \cos (c+d x)+a}}\right )}{\sqrt {a} d}+\frac {2 (31 A-7 B+35 C) \sin (c+d x) \sec ^{\frac {3}{2}}(c+d x)}{105 d \sqrt {a \cos (c+d x)+a}}-\frac {2 (43 A-91 B+35 C) \sin (c+d x) \sqrt {\sec (c+d x)}}{105 d \sqrt {a \cos (c+d x)+a}}-\frac {2 (A-7 B) \sin (c+d x) \sec ^{\frac {5}{2}}(c+d x)}{35 d \sqrt {a \cos (c+d x)+a}}+\frac {2 A \sin (c+d x) \sec ^{\frac {7}{2}}(c+d x)}{7 d \sqrt {a \cos (c+d x)+a}} \]

[In]

Int[((A + B*Cos[c + d*x] + C*Cos[c + d*x]^2)*Sec[c + d*x]^(9/2))/Sqrt[a + a*Cos[c + d*x]],x]

[Out]

(Sqrt[2]*(A - B + C)*ArcTan[(Sqrt[a]*Sin[c + d*x])/(Sqrt[2]*Sqrt[Cos[c + d*x]]*Sqrt[a + a*Cos[c + d*x]])]*Sqrt
[Cos[c + d*x]]*Sqrt[Sec[c + d*x]])/(Sqrt[a]*d) - (2*(43*A - 91*B + 35*C)*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(105
*d*Sqrt[a + a*Cos[c + d*x]]) + (2*(31*A - 7*B + 35*C)*Sec[c + d*x]^(3/2)*Sin[c + d*x])/(105*d*Sqrt[a + a*Cos[c
 + d*x]]) - (2*(A - 7*B)*Sec[c + d*x]^(5/2)*Sin[c + d*x])/(35*d*Sqrt[a + a*Cos[c + d*x]]) + (2*A*Sec[c + d*x]^
(7/2)*Sin[c + d*x])/(7*d*Sqrt[a + a*Cos[c + d*x]])

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 2861

Int[1/(Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*Sqrt[(c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> D
ist[-2*(a/f), Subst[Int[1/(2*b^2 - (a*c - b*d)*x^2), x], x, b*(Cos[e + f*x]/(Sqrt[a + b*Sin[e + f*x]]*Sqrt[c +
 d*Sin[e + f*x]]))], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 -
 d^2, 0]

Rule 3063

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[(B*c - A*d)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x]
)^(n + 1)/(f*(n + 1)*(c^2 - d^2))), x] + Dist[1/(b*(n + 1)*(c^2 - d^2)), Int[(a + b*Sin[e + f*x])^m*(c + d*Sin
[e + f*x])^(n + 1)*Simp[A*(a*d*m + b*c*(n + 1)) - B*(a*c*m + b*d*(n + 1)) + b*(B*c - A*d)*(m + n + 2)*Sin[e +
f*x], x], x], x] /; FreeQ[{a, b, c, d, e, f, A, B, m}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2
 - d^2, 0] && LtQ[n, -1] && (IntegerQ[n] || EqQ[m + 1/2, 0])

Rule 3122

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_)*((A_.) + (B_.)*s
in[(e_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> Simp[(-(c^2*C - B*c*d + A*d^2))*Cos[e
+ f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^(n + 1)/(d*f*(n + 1)*(c^2 - d^2))), x] + Dist[1/(b*d*(n +
1)*(c^2 - d^2)), Int[(a + b*Sin[e + f*x])^m*(c + d*Sin[e + f*x])^(n + 1)*Simp[A*d*(a*d*m + b*c*(n + 1)) + (c*C
 - B*d)*(a*c*m + b*d*(n + 1)) + b*(d*(B*c - A*d)*(m + n + 2) - C*(c^2*(m + 1) + d^2*(n + 1)))*Sin[e + f*x], x]
, x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C, m}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^
2, 0] &&  !LtQ[m, -2^(-1)] && (LtQ[n, -1] || EqQ[m + n + 2, 0])

Rule 4306

Int[(u_)*((c_.)*sec[(a_.) + (b_.)*(x_)])^(m_.), x_Symbol] :> Dist[(c*Sec[a + b*x])^m*(c*Cos[a + b*x])^m, Int[A
ctivateTrig[u]/(c*Cos[a + b*x])^m, x], x] /; FreeQ[{a, b, c, m}, x] &&  !IntegerQ[m] && KnownSineIntegrandQ[u,
 x]

Rubi steps \begin{align*} \text {integral}& = \left (\sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx \\ & = \frac {2 A \sec ^{\frac {7}{2}}(c+d x) \sin (c+d x)}{7 d \sqrt {a+a \cos (c+d x)}}+\frac {\left (2 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {-\frac {1}{2} a (A-7 B)+\frac {1}{2} a (6 A+7 C) \cos (c+d x)}{\cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx}{7 a} \\ & = -\frac {2 (A-7 B) \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d \sqrt {a+a \cos (c+d x)}}+\frac {2 A \sec ^{\frac {7}{2}}(c+d x) \sin (c+d x)}{7 d \sqrt {a+a \cos (c+d x)}}+\frac {\left (4 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {\frac {1}{4} a^2 (31 A-7 B+35 C)-a^2 (A-7 B) \cos (c+d x)}{\cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx}{35 a^2} \\ & = \frac {2 (31 A-7 B+35 C) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d \sqrt {a+a \cos (c+d x)}}+\frac {2 A \sec ^{\frac {7}{2}}(c+d x) \sin (c+d x)}{7 d \sqrt {a+a \cos (c+d x)}}+\frac {\left (8 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {-\frac {1}{8} a^3 (43 A-91 B+35 C)+\frac {1}{4} a^3 (31 A-7 B+35 C) \cos (c+d x)}{\cos ^{\frac {3}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx}{105 a^3} \\ & = -\frac {2 (43 A-91 B+35 C) \sqrt {\sec (c+d x)} \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d \sqrt {a+a \cos (c+d x)}}+\frac {2 A \sec ^{\frac {7}{2}}(c+d x) \sin (c+d x)}{7 d \sqrt {a+a \cos (c+d x)}}+\frac {\left (16 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {105 a^4 (A-B+C)}{16 \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}} \, dx}{105 a^4} \\ & = -\frac {2 (43 A-91 B+35 C) \sqrt {\sec (c+d x)} \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d \sqrt {a+a \cos (c+d x)}}+\frac {2 A \sec ^{\frac {7}{2}}(c+d x) \sin (c+d x)}{7 d \sqrt {a+a \cos (c+d x)}}+\left ((A-B+C) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {1}{\sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}} \, dx \\ & = -\frac {2 (43 A-91 B+35 C) \sqrt {\sec (c+d x)} \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d \sqrt {a+a \cos (c+d x)}}+\frac {2 A \sec ^{\frac {7}{2}}(c+d x) \sin (c+d x)}{7 d \sqrt {a+a \cos (c+d x)}}-\frac {\left (2 a (A-B+C) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \text {Subst}\left (\int \frac {1}{2 a^2+a x^2} \, dx,x,-\frac {a \sin (c+d x)}{\sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}\right )}{d} \\ & = \frac {\sqrt {2} (A-B+C) \arctan \left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}}{\sqrt {a} d}-\frac {2 (43 A-91 B+35 C) \sqrt {\sec (c+d x)} \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{105 d \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d \sqrt {a+a \cos (c+d x)}}+\frac {2 A \sec ^{\frac {7}{2}}(c+d x) \sin (c+d x)}{7 d \sqrt {a+a \cos (c+d x)}} \\ \end{align*}

Mathematica [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 5 vs. order 3 in optimal.

Time = 8.21 (sec) , antiderivative size = 2607, normalized size of antiderivative = 10.14 \[ \int \frac {\left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\text {Result too large to show} \]

[In]

Integrate[((A + B*Cos[c + d*x] + C*Cos[c + d*x]^2)*Sec[c + d*x]^(9/2))/Sqrt[a + a*Cos[c + d*x]],x]

[Out]

(2*Cos[c/2 + (d*x)/2]*Sqrt[(1 - 2*Sin[c/2 + (d*x)/2]^2)^(-1)]*Sqrt[1 - 2*Sin[c/2 + (d*x)/2]^2]*((2*B*Sin[c/2 +
 (d*x)/2])/(7*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(7/2)) - (2*C*Sin[c/2 + (d*x)/2])/(7*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(
7/2)) + (12*B*Sin[c/2 + (d*x)/2])/(35*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(5/2)) + (2*C*Sin[c/2 + (d*x)/2])/(35*(1 -
2*Sin[c/2 + (d*x)/2]^2)^(5/2)) + (16*B*Sin[c/2 + (d*x)/2])/(35*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(3/2)) + (8*C*Sin[
c/2 + (d*x)/2])/(105*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(3/2)) + (32*B*Sin[c/2 + (d*x)/2])/(35*Sqrt[1 - 2*Sin[c/2 +
(d*x)/2]^2]) + (16*C*Sin[c/2 + (d*x)/2])/(105*Sqrt[1 - 2*Sin[c/2 + (d*x)/2]^2]) + ((A - B + C)*Csc[c/2 + (d*x)
/2]^9*(363825*Sin[c/2 + (d*x)/2]^2 - 4729725*Sin[c/2 + (d*x)/2]^4 + 26785605*Sin[c/2 + (d*x)/2]^6 - 86790165*S
in[c/2 + (d*x)/2]^8 + 177677808*Sin[c/2 + (d*x)/2]^10 - 239283044*Sin[c/2 + (d*x)/2]^12 + 52080*Hypergeometric
2F1[2, 11/2, 13/2, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^12 + 560*Hypergeomet
ricPFQ[{2, 2, 2, 2, 11/2}, {1, 1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x
)/2]^12 + 213120160*Sin[c/2 + (d*x)/2]^14 - 168280*Hypergeometric2F1[2, 11/2, 13/2, Sin[c/2 + (d*x)/2]^2/(-1 +
 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^14 - 2240*HypergeometricPFQ[{2, 2, 2, 2, 11/2}, {1, 1, 1, 13/2},
Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^14 - 121497024*Sin[c/2 + (d*x)/2]^16 +
212520*Hypergeometric2F1[2, 11/2, 13/2, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]
^16 + 3360*HypergeometricPFQ[{2, 2, 2, 2, 11/2}, {1, 1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)
/2]^2)]*Sin[c/2 + (d*x)/2]^16 + 40125184*Sin[c/2 + (d*x)/2]^18 - 124320*Hypergeometric2F1[2, 11/2, 13/2, Sin[c
/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^18 - 2240*HypergeometricPFQ[{2, 2, 2, 2, 11/
2}, {1, 1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^18 - 5840384*Sin[c
/2 + (d*x)/2]^20 + 28000*Hypergeometric2F1[2, 11/2, 13/2, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*
Sin[c/2 + (d*x)/2]^20 + 560*HypergeometricPFQ[{2, 2, 2, 2, 11/2}, {1, 1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 +
2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^20 + 363825*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d
*x)/2]^2)]]*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] - 5336100*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]
^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^2*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2
)] + 34636140*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^4*Sqrt[Sin[
c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] - 131060160*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2
+ (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^6*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] + 320535600*ArcT
anh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^8*Sqrt[Sin[c/2 + (d*x)/2]^2/(
-1 + 2*Sin[c/2 + (d*x)/2]^2)] - 530671680*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Si
n[c/2 + (d*x)/2]^10*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] + 604296000*ArcTanh[Sqrt[Sin[c/2
+ (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^12*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 +
 (d*x)/2]^2)] - 468948480*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]
^14*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] + 237726720*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1
 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^16*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] -
70963200*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^18*Sqrt[Sin[c/2
+ (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] + 9461760*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x
)/2]^2)]]*Sin[c/2 + (d*x)/2]^20*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] - 1120*Cos[(c + d*x)/
2]^6*HypergeometricPFQ[{2, 2, 2, 11/2}, {1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[
c/2 + (d*x)/2]^12*(-6 + 5*Sin[c/2 + (d*x)/2]^2) + 280*Cos[(c + d*x)/2]^4*HypergeometricPFQ[{2, 2, 11/2}, {1, 1
3/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^12*(103 - 164*Sin[c/2 + (d*x)/2]^
2 + 70*Sin[c/2 + (d*x)/2]^4)))/(40425*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(9/2)*(-1 + 2*Sin[c/2 + (d*x)/2]^2))))/(d*S
qrt[a*(1 + Cos[c + d*x])])

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(496\) vs. \(2(218)=436\).

Time = 1.86 (sec) , antiderivative size = 497, normalized size of antiderivative = 1.93

method result size
default \(-\frac {\sqrt {2}\, \left (\sec ^{\frac {9}{2}}\left (d x +c \right )\right ) \sqrt {\left (1+\cos \left (d x +c \right )\right ) a}\, \left (105 A \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{5}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}-105 \left (\cos ^{5}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, B \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )+105 C \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{5}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}+105 A \left (\cos ^{4}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )+43 A \left (\cos ^{4}\left (d x +c \right )\right ) \sin \left (d x +c \right ) \sqrt {2}-105 B \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{4}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}-91 \left (\cos ^{4}\left (d x +c \right )\right ) \sin \left (d x +c \right ) B \sqrt {2}+105 C \left (\cos ^{4}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )+35 C \left (\cos ^{4}\left (d x +c \right )\right ) \sin \left (d x +c \right ) \sqrt {2}-31 A \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right ) \sqrt {2}+7 \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right ) B \sqrt {2}-35 C \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right ) \sqrt {2}+3 A \left (\cos ^{2}\left (d x +c \right )\right ) \sin \left (d x +c \right ) \sqrt {2}-21 \left (\cos ^{2}\left (d x +c \right )\right ) \sin \left (d x +c \right ) B \sqrt {2}-15 A \cos \left (d x +c \right ) \sin \left (d x +c \right ) \sqrt {2}\right )}{105 d a \left (1+\cos \left (d x +c \right )\right )}\) \(497\)
parts \(-\frac {A \sqrt {2}\, \left (\sec ^{\frac {9}{2}}\left (d x +c \right )\right ) \sqrt {\left (1+\cos \left (d x +c \right )\right ) a}\, \left (105 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{5}\left (d x +c \right )\right )+105 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{4}\left (d x +c \right )\right )+43 \sqrt {2}\, \left (\cos ^{4}\left (d x +c \right )\right ) \sin \left (d x +c \right )-31 \sqrt {2}\, \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right )+3 \sqrt {2}\, \left (\cos ^{2}\left (d x +c \right )\right ) \sin \left (d x +c \right )-15 \sqrt {2}\, \cos \left (d x +c \right ) \sin \left (d x +c \right )\right )}{105 d a \left (1+\cos \left (d x +c \right )\right )}+\frac {B \sqrt {2}\, \left (\sec ^{\frac {9}{2}}\left (d x +c \right )\right ) \sqrt {\left (1+\cos \left (d x +c \right )\right ) a}\, \left (15 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{5}\left (d x +c \right )\right )+15 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{4}\left (d x +c \right )\right )+13 \sqrt {2}\, \left (\cos ^{4}\left (d x +c \right )\right ) \sin \left (d x +c \right )-\sqrt {2}\, \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right )+3 \sqrt {2}\, \left (\cos ^{2}\left (d x +c \right )\right ) \sin \left (d x +c \right )\right )}{15 d a \left (1+\cos \left (d x +c \right )\right )}-\frac {C \sqrt {2}\, \left (\sec ^{\frac {9}{2}}\left (d x +c \right )\right ) \sqrt {\left (1+\cos \left (d x +c \right )\right ) a}\, \left (3 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{5}\left (d x +c \right )\right )+3 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{4}\left (d x +c \right )\right )+\sqrt {2}\, \left (\cos ^{4}\left (d x +c \right )\right ) \sin \left (d x +c \right )-\sqrt {2}\, \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right )\right )}{3 d a \left (1+\cos \left (d x +c \right )\right )}\) \(569\)

[In]

int((A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^(9/2)/(a+a*cos(d*x+c))^(1/2),x,method=_RETURNVERBOSE)

[Out]

-1/105/d*2^(1/2)/a*sec(d*x+c)^(9/2)*((1+cos(d*x+c))*a)^(1/2)/(1+cos(d*x+c))*(105*A*arcsin(cot(d*x+c)-csc(d*x+c
))*cos(d*x+c)^5*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)-105*cos(d*x+c)^5*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)*B*arcsin(
cot(d*x+c)-csc(d*x+c))+105*C*arcsin(cot(d*x+c)-csc(d*x+c))*cos(d*x+c)^5*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)+105*
A*cos(d*x+c)^4*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)*arcsin(cot(d*x+c)-csc(d*x+c))+43*A*cos(d*x+c)^4*sin(d*x+c)*2^
(1/2)-105*B*arcsin(cot(d*x+c)-csc(d*x+c))*cos(d*x+c)^4*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)-91*cos(d*x+c)^4*sin(d
*x+c)*B*2^(1/2)+105*C*cos(d*x+c)^4*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)*arcsin(cot(d*x+c)-csc(d*x+c))+35*C*cos(d*
x+c)^4*sin(d*x+c)*2^(1/2)-31*A*cos(d*x+c)^3*sin(d*x+c)*2^(1/2)+7*cos(d*x+c)^3*sin(d*x+c)*B*2^(1/2)-35*C*cos(d*
x+c)^3*sin(d*x+c)*2^(1/2)+3*A*cos(d*x+c)^2*sin(d*x+c)*2^(1/2)-21*cos(d*x+c)^2*sin(d*x+c)*B*2^(1/2)-15*A*cos(d*
x+c)*sin(d*x+c)*2^(1/2))

Fricas [A] (verification not implemented)

none

Time = 0.31 (sec) , antiderivative size = 189, normalized size of antiderivative = 0.74 \[ \int \frac {\left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=-\frac {\frac {105 \, \sqrt {2} {\left ({\left (A - B + C\right )} a \cos \left (d x + c\right )^{4} + {\left (A - B + C\right )} a \cos \left (d x + c\right )^{3}\right )} \arctan \left (\frac {\sqrt {2} \sqrt {a \cos \left (d x + c\right ) + a} \sqrt {\cos \left (d x + c\right )}}{\sqrt {a} \sin \left (d x + c\right )}\right )}{\sqrt {a}} + \frac {2 \, {\left ({\left (43 \, A - 91 \, B + 35 \, C\right )} \cos \left (d x + c\right )^{3} - {\left (31 \, A - 7 \, B + 35 \, C\right )} \cos \left (d x + c\right )^{2} + 3 \, {\left (A - 7 \, B\right )} \cos \left (d x + c\right ) - 15 \, A\right )} \sqrt {a \cos \left (d x + c\right ) + a} \sin \left (d x + c\right )}{\sqrt {\cos \left (d x + c\right )}}}{105 \, {\left (a d \cos \left (d x + c\right )^{4} + a d \cos \left (d x + c\right )^{3}\right )}} \]

[In]

integrate((A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^(9/2)/(a+a*cos(d*x+c))^(1/2),x, algorithm="fricas")

[Out]

-1/105*(105*sqrt(2)*((A - B + C)*a*cos(d*x + c)^4 + (A - B + C)*a*cos(d*x + c)^3)*arctan(sqrt(2)*sqrt(a*cos(d*
x + c) + a)*sqrt(cos(d*x + c))/(sqrt(a)*sin(d*x + c)))/sqrt(a) + 2*((43*A - 91*B + 35*C)*cos(d*x + c)^3 - (31*
A - 7*B + 35*C)*cos(d*x + c)^2 + 3*(A - 7*B)*cos(d*x + c) - 15*A)*sqrt(a*cos(d*x + c) + a)*sin(d*x + c)/sqrt(c
os(d*x + c)))/(a*d*cos(d*x + c)^4 + a*d*cos(d*x + c)^3)

Sympy [F(-1)]

Timed out. \[ \int \frac {\left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\text {Timed out} \]

[In]

integrate((A+B*cos(d*x+c)+C*cos(d*x+c)**2)*sec(d*x+c)**(9/2)/(a+a*cos(d*x+c))**(1/2),x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {\left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\int { \frac {{\left (C \cos \left (d x + c\right )^{2} + B \cos \left (d x + c\right ) + A\right )} \sec \left (d x + c\right )^{\frac {9}{2}}}{\sqrt {a \cos \left (d x + c\right ) + a}} \,d x } \]

[In]

integrate((A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^(9/2)/(a+a*cos(d*x+c))^(1/2),x, algorithm="maxima")

[Out]

integrate((C*cos(d*x + c)^2 + B*cos(d*x + c) + A)*sec(d*x + c)^(9/2)/sqrt(a*cos(d*x + c) + a), x)

Giac [F]

\[ \int \frac {\left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\int { \frac {{\left (C \cos \left (d x + c\right )^{2} + B \cos \left (d x + c\right ) + A\right )} \sec \left (d x + c\right )^{\frac {9}{2}}}{\sqrt {a \cos \left (d x + c\right ) + a}} \,d x } \]

[In]

integrate((A+B*cos(d*x+c)+C*cos(d*x+c)^2)*sec(d*x+c)^(9/2)/(a+a*cos(d*x+c))^(1/2),x, algorithm="giac")

[Out]

integrate((C*cos(d*x + c)^2 + B*cos(d*x + c) + A)*sec(d*x + c)^(9/2)/sqrt(a*cos(d*x + c) + a), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {\left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right ) \sec ^{\frac {9}{2}}(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\int \frac {{\left (\frac {1}{\cos \left (c+d\,x\right )}\right )}^{9/2}\,\left (C\,{\cos \left (c+d\,x\right )}^2+B\,\cos \left (c+d\,x\right )+A\right )}{\sqrt {a+a\,\cos \left (c+d\,x\right )}} \,d x \]

[In]

int(((1/cos(c + d*x))^(9/2)*(A + B*cos(c + d*x) + C*cos(c + d*x)^2))/(a + a*cos(c + d*x))^(1/2),x)

[Out]

int(((1/cos(c + d*x))^(9/2)*(A + B*cos(c + d*x) + C*cos(c + d*x)^2))/(a + a*cos(c + d*x))^(1/2), x)